Intermittency in two dimensions 
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We introduce a family of area-preserving maps representing a (non-trivial) two-dimensional ex- 
tension of the Pomeau-Manneville family in one dimension. We analyze the long-time behavior of 
recurrence time distributions and correlations, providing analytical and numerical estimates. We 
study the transport properties of a suitable lift and use a probabilistic argument to derive the full 
spectrum of transport moments. Finally the dynamical effects of a stochastic perturbation are 
considered. 
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I. INTRODUCTION 

Hamiltonian systems arc gcnerically not fully hyper- 
bolic [1]: for example the phase space of typical area- 
preserving maps reveals the co-existence of chaotic tra- 
jectories and islands of regular motion (periodic or quasi- 
periodic trajectories) 0]. Even when we are concerned 
with statistical properties of motion on the chaotic com- 
ponent we cannot neglect the presence of regular struc- 
tures. They deeply influence chaotic motion as, when- 
ever trajectories come close to integrable islands, they 
stick there for some time and irregular dynamics is thus 
punctuated by laminar segments where the system 'mim- 
ics' an integrable one. This intermittent behavior has 
strong influence on the long-time properties of quanti- 
ties like correlations decay or recurrence time statistics, 
that typically present a power-law tail [H 0, [H H, 0|- 
In unbounded systems intermittency influences transport 
properties, generating anomalous diffusion processes (see 
8] and references therein) , in contrast to the normal dif- 
fusion observed for fully hyperbolic systems 0]. 

While much effort is still devoted to fully understand 
the general picture of mixed phase space, the situation 
simplifies if we let the islands of regular motion shrink 
to zero: even the presence of a single marginally sta- 
ble fixed point can produce intermittent-like behavior 
[lOl . Ull . In one dimension this corresponds to the 

Pomeau-Manneville maps on the unit interval [l3| 



°n\mod 1 



Z > 1 



(1) 



which represent one of the few examples of non-fully hy- 
perbolic systems for which analytic results can be ob- 
tained with a variety of different techniques (see for in- 
stance [HI, EH EH E3, EM EH EH Hl| ) - Such maps present 
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a polynomial decay of correlations and recurrence time 
statistics with exponents that depend on the intermit- 
tency parameter 7 [TEl EH [S3- Moreover a proper lift 
on the real line can generate anomalous diffusion and the 
set of transport moments typically shows a two-scales 
structure [Si El EH El] ■ 

The situation is much less satisfactory in more than 
one dimension. Rigorous results were derived for spe- 
cific cases, where it was possible to give precise bounds 
on the rate of mixing [lj, Ell- Here we study a fam- 
ily of area- preserving maps with a neutral fixed point. 
This family depends on a parameter that governs sta- 
bility properties of the fixed point, in analogy with the 
Pomeau-Manneville maps. 

We introduce the two dimensional family of area- 
preserving maps in section II, where we also discuss the 
dynamical features we will look at. In section III wc 
will consider the unstable manifold of the neutral fixed 
point, and provide simple estimates that will be pivotal 
in predicting the decay of survival probabilities. Section 
IV contains extensive investigations on survival probabil- 
ity and correlation functions decay. By lifting the map 
on an unbounded phase space we then analyze transport 
properties in section V. The role of a stochastic pertur- 
bation is then studied in section VI, while we present our 
conclusions in section VII. 



II. THE MODEL 

We define the following one-parameter family of maps 
T 7 (x,y) : T 2 -> T 2 , where T 2 = [-7r,7r) 2 (with torus 
topology): 



T 7 (x,y) 



x + f-y(x) + y on T 
y + f-y (x) on T 



with 



/ 7 (x) = 7rsign(a;) 



7 > 1. 



(2) 



(3) 



2 



The map T 7 is area-preserving for every choice of the 
impulsive force f(x). The Jacobian matrix is 



l + 1 

m i 



(4) 



so we have det J 7 (x, y) = 1 and Tr J 7 (x, y) = 2 + /'(x) > 
2 for x G [— 7r, 7r)/{0}: the map is everywhere hyperbolic 
except on the line x — and thus the fixed point at 
the origin is marginally stable (parabolic) for every 7. 
The parameter 7 changes the 'stickiness' of the origin, in 
analogy with the intermittency exponent z appearing in 
the Pomeau-Manneville maps of E g. C D). We point out 
that the choice / (x) — x — sin(x) [l(J, LL2| gives rise to 
a marginal fixed point with the same stability properties 
of that for 7 = 3 in Eq. £2]). We note explicitly that for 
our model / € C k with k — [7 — 1] (where [•] denotes the 
integer part) unless 7 is an odd integer, in which case 

/ G c°° H]. 



A. Dynamical indicators 

In order to get information about the dynamical prop- 
erties of the systems it is often useful to employ time 
statistics [2!| . We choose a set f2 including the parabolic 
fixed point and then define a partition of fi in disjoint 
sets fl n , each representing the set of points that leave Q 
in exactly n iterations. The survival probability pn(n) 
is the fraction of initial conditions in fl that are still in 
f2 after n iterations. The behavior of pn(n) generally 
depends on the choice of the measure \ii with which we 
distribute initial conditions over f2, which may be quite 
different from the invariant measure. In the present case 
the invariant measure is the Lebesgue one fj,, which also 
represents the most natural way to spread initial condi- 
tions over Q, so /ii — fj, and 



Pn(n) 



1 



(5) 



We may also define the waiting time distribution (or res- 
idence time statistics) over fi, ^n(n), as the probability 
that once a trajectory enters the set fl it stays there ex- 
actly n time steps. ipn(n) is computed by running a long 
trajectory and recording residence times in Q: ip(i( n ) is 
just the probability distribution of such residence times. 
In our case 



I 



(/x(fl n ) - ji(O n+ i)) 



(6) 



where ergodicity and the property that the map preserves 
Lebesgue measure have been used. We point out that in 
general, while pn(n) depends upon an -arbitrary- choice 
of the distribution of initial conditions, il>a(n) doesn't. 



From Eq. C51 [6]) we see that the asymptotics of these 
quantities are tightly related (29j : in particular if the 
measure of the sets fl n decays according to a power law 



we get 



and 



Pn{n) 



(7) 



(8) 



(9) 



Besides their intrinsic interest, these relations bear re- 
markable links with the problem of establishing the mix- 
ing rates for the system 0, M, US HH • We briefly recall 
it with a simple argument: suppose we consider an ob- 
servable A that remains fully correlated for portions of 
trajectories within f2, and otherwise completely decorre- 
lated (due to randomness of motion outside of CI); then 
we may easily show that [3l[ 

C AA (n) = (A(n)A(0)) - (A) 2 (10) 

~ ((A 2 ) - (A) 2 ) £ jr Mk); 



m=n k=m 



so that the exponents of power-law decay of survival 
probability and correlations should coincide 



C AA {n) ~ n- a . 



(11) 



We observe that the validity of Eq. (fTU|) has been care- 
fully numerically scrutinized for chaotic billiards [3l|, HH , 
and even validated in rigorous estimates of polynomial 
mixing speed for Id intermittent systems [Iff (see also 
[33J), but indications of its possible failure have also been 
suggested [34| . 

It is interesting to remark that showing that power law 
decays of pn(n) and ipn(n) differ by two (Eq. (JS1 E])) em- 
ploys the fact that Lebesgue measure is the invariant one 
for the system. Actually for the map of Eq. £!]) , where 
the invariant measure is not uniform, they differ by one 
if we choose the initial conditions uniformly distributed 
with Lebesgue measure (while the exponent of the wait- 
ing time distribution coincides with the one ruling the 
length of the segments fi ra [IH) . A relationship coincid- 
ing with Eq. (0 [5]) holds instead for another intermittent 
map, introduced by Pikovsky [2~j ] (some features of this 
map are also described in [35j] ) : 



1+1 = fz{x n ), 



(12) 



where f z is an (odd) circle map, again dependent on 
an intermittency parameter z, implicitly defined on the 
torus T = [—1, 1) by 



/*(*) + e(i 



< x < 1/(22) 
l/(2z) < x < 1 



(13) 
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A key feature that the map of Eq. (fT2| shares with our 
model is that the invariant distribution is smooth, coin- 
ciding with the Lebesgue measure, as it can be checked 
by inspection of the form of Perron- Frobenius operator. 

For the map under consideration it is possible to get 
an estimate of the exponent a by studying the invariant 
manifolds of the parabolic fixed point. 



III. INVARIANT MANIFOLDS 

A typical trajectory staying for a long time in f2 (again 
we take as Q a region including the parabolic fixed point) 
enters f2 close to the stable manifold, escaping along the 
unstable one (Fig. (Q])). In particular we are going to 
discuss how trajectories escape by following the unstable 
manifold of the marginal fixed point. For the odd sym- 
metry of fj(x), we can restrict the analysis to the first 
quadrant. 



0.25 




-0.25 



FIG. 1: (color online) Portion of phase space of the map of Eq. 
© for 7 = 3 close to the marginal fixed point together with 
the graph of its unstable manifold (continuous (red) line). 

Let's call (x,y(x)) the graph of the unstable manifold 
in the neighborhood of the origin and suppose that, very 
close to the indifferent fixed point, y(x) ~ x° . With the 
choice of Eq. we have that f 1 (x) ~ a x 1 and the 
Jacobian matrix (Eq. of the map is: 



J 7 (x,y) 



1 + b x^- 1 
bxi' 1 



(14) 



whose eigenvalues are written to leading order as: X± = 
l±Px^- 1 ^ 2 . 

The vector (l,y'(x)), i.e. (l,r) x a ~ 1 ), tangent to the 
unstable manifold satisfies 



J 7 (x,y) 



(^-)" (1 + /3 ^ 1)/2) (^ Ll 



(15) 



1 + 6 x~i~ x + r] 



bx^- 1 



r\ x 



~ rjx*- 1 + -q!3 x { ' 1 - 1 ^ 2+{ ' 7 -^ 

and from these equations, remembering that x « 1 and 
7 > 1, we obtain 



7+1 



(17) 



We note explicitely that for 7 = 3, this result is in agree- 
ment with the case / = x — sin(x) derived in [jjj]. We 
can now study the dynamics restricted to the unstable 
manifold. Let's call £ the arclength coordinate along the 
manifold; for small x we get: 



£{x) 



dxsjl + (dy(x)/dx) 2 ~ x (18) 



Denote by £ n the coordinate £ at a point (x n ,y(x n )) 
and by £ n +i the coordinate along the manifold of 

X!y {x n: y(x n )) 



4+1 = in + h(£ n ). 
By using Eq. (JTBJ) and © we get 



(19) 



Hi) 



d( 
~dt 



d£ dx 
dx dt 



(£) ~ {y{£)+x^(£)) =r+P (20) 



From Eq. I|17p and by remembering that 7 > 1 we obtain, 
via a continuous time approximation (36l . [37| , 



h(£) 



dt' 



(21) 



If we fix the boundary of £1 at a scale L we can then eval- 
uate the time needed to reach the boundary as a function 
of the arclength £ along the manifold, by employing the 
standard argument of [3q . \3l\ : 



T(£) 



(22) 



that implies the following scaling for the inverse function 
£{T) - T~^+ (23) 

We now arrive to the crucial point: we estimate pn(n) 
as the area of rectangle having one vertex at the origin 
(the parabolic fixed point), and another at a point on the 
unstable manifold (x, y) that exits f2 in n steps (see Fig. 
©): 



x-y 



£{n) 



<T + 1 



(n~ — Y +1 = n~^r (24) 



Thus we have an estimate of power law decay of the 
survival probability as a function of the intermittency 
parameter 7 as 



Pn(n) 



-1+3 
n -t- 1 



(25) 
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Po( n ) 




FIG. 2: (color online) A few Q n (once we set fi as the first 
quadrant x > 0, y > 0). 



and for the waiting time distribution as well (Eq. 

37 + 1 

Wn{n) ~ n . 



(26) 



In view of the argument we earlier mentioned (see Eq. 
JTTJ)), the estimate of Eq. ([23]) suggests the same decay 
law for (auto)correlation functions 



(27) 



Next section will present several numerical simulations 
concerning these quantities. 



IV. ASYMPTOTIC DECAYS 

We start by considering the survival probability: Fig. 
(|3|) shows two examples of numerically computed Pfi(n). 
The numerical data exhibit an excellent agreement with 
analytic estimates over a wide range of intermittency pa- 
rameters, as shown in Fig. ((4]), which also provides clear 
indications of the validity of our estimate for the asymp- 
totic decay of the waiting time distribution. 

We already mentioned that various arguments support 
the expectation that correlation functions should decay 
as the survival probability (Eq. (|27|) ). so we proceed 
to scrutinize this prediction by running extensive direct 
numerical simulations on autocorrelation functions; as 
we are dealing with an ergodic (and mixing pjj]) system, 
autocorrelation functions can be evaluated in terms of 
phase space averages (instead of temporal averages): 



CAA(n) 



d[i(z)A(T"z)A(z) - 



M 



dfi(z)A(z) 



M 



FIG. 3: (color online) Survival probabilities for 7 = 3 (lower 
curve) and 7 = 10 (upper curve) together with the power law 
decays predicted by Eq. (|25p . We used 10 12 initial conditions 
and set Q = [- 



1 in 2 
2 ' 21 



12 r 




(28) 



FIG. 4: (color online) Exponents of power-law decay for sur- 
vival probability pn{n) and waiting time distribution ipn(n): 
lines refer to analytic estimates of Eq. (|25[) -upper. (126[) -lower: 
open circles come from numerical simulations of the waiting 
time distribution, open diamonds from survival probability 
simulations. 



where A is a smooth function on the phase space M. 
and n is the invariant Lebesgue measure. From a nu- 
merical point of view it is known that often Monte-Carlo 
evaluation of Eq. (f28j) cannot be pushed too far, as the 
statistical error is of order 1/s/N in the number of ini- 
tial conditions: generally we also expect an (exponential) 
transient in the decay [13, H3, 38] : transient time t might 
depend on both 7 and the choice of phase space function 
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A [6[ . We also remark that smoothness of the function 
A plays a fundamental role: as a matter of fact we may 
obtain arbitrarily slow correlation decay even for Anosov 
maps by using integrable non-smooth functions [3^ |. or, 
from a mathematical point of view, we may have that 
the degree of smoothness determines the essential spec- 
tral radius of the Perron- Frobenius operator ■ 

We performed the explicit calculation of the autocor- 
relation function for different values of the intermittency 
parameter 7 and for different observables. We obtained 
the best results (i.e. cleanest curves and shortest time t) 
for large values of 7 and by using A(x, y) = e~ v . The 
choice of the smooth function to use is quite arbitrary; 
we looked for a function not vanishing in correspondence 
of the marginal fixed point (as suggested for example in 
[H) and by the special choice of a gaussian depending on 
a single variable we could save computational time (see 
also |2l|). In figures (JSJ and (O we present results for 
7 = 3 and 7 = 10: 




FIG. 5: (color online) Autocorrelation function for the ob- 
servable e~ v and 7 = 3. We used 2 • 10 10 initial conditions 
(uniformly distributed in the torus cell). The predicted decay 
is shown by the (red) straight line, which has a slope —3. 

The case reported in Fig. ([5]) is important because 
for a similar 2d mapping (having the same intermittent 
exponent 7 = 3) it was proved in [12j that the decay 
is faster than n~ 2 , and a class of cross-correlation was 
constructed indicating that the bound is optimal: while 
survival probability data for the same exponent indicate 
clearly that the decay we predict (n~ 3 ) is numerically 
well reproduced, data for correlations are less conclusive 
(see Fig. ([5])). In general, numerical data are more dif- 
ficult to interpret for low values of 7, and numerical fits 
tend to lie below the predicted exponents (see Fig. |J7J), 
while for larger values of 7 the accordance with our esti- 
mates is much better. Moreover, the agreement improves 
by increasing the number of initial conditions, so that we 




FIG. 6: (color online) Autocorrelation function for the ob- 

— 2 Q 
servable e y and 7 = 10. We used 10 initial conditions 

(uniformly distributed in the torus cell). The predicted decay 

is shown by the (red) straight line, which has a slope —13/9. 



expect the discrepancy to be essentially due to numerical 
limitations (Fig. ||7J)). Further indications of the similar- 




FIG. 7: (color online) Numerical values of power-law decay 

exponents for the autocorrelation function for the observable 
_ 2 

e v ((black) circles and (blue) triangles) together with the 
analytical estimates (full (red) line). Circles were obtained by 
using 2-10 9 initial conditions while the triangles were obtained 
for 7 = 2.5 by using 5 ■ 10 10 initial conditions, in the range 
3 < 7 < 6 by using 2 • 10 10 initial conditions, and in the range 
6.5 < 7 < 10 by using 10 10 initial conditions. 

ity between correlations and survival probability decays 
will be provided in section VI, when considering the role 
of stochastic perturbations. 
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V. TRANSPORT PROPERTIES 

In order to study transport properties, we have to 
abandon the dynamics restricted to the torus (Eq. @) 
by lifting the map in an appropriate way. 

For the sake of clarity we introduce a third dimension 
(say z) to describe the motion through elementary cells. 
We then assign a jumping number +1 to the points be- 
longing to the first quadrant, —I to the points belonging 
to the third quadrant and to all the other points. This 
means that a laminar phase of length n will correspond 
to a jump in the positive direction of n elementary cells 
(Fig®. 



+3T 



-Jt 






+1 


-1 






m+1 




m-1 



FIG. 8: Jumping numbers and lifted map. 
Formally, the lift is given by the following formula: 



T 7 (x,y,m) = 



(T 7 {x,y),m) 



for xy < 



(T^(x, y), m + sign(x)) forxy > 



,(29) 



where m is an integer variable. 

Considering successive entrance in the laminar regions 
as uncorrelated we can approximate the diffusion process 
by a Continuous Time Random Walk (CTRW) [23, 
with the probability distribution of the laminar phases 
given by the waiting time distribution ipa{n). 
In particular, by making use of the CTRW approach it 
is possible to characterize the transport properties of the 
process in terms of the set of moments of the diffusing 
variable [43| : 

(|z(n)-z(0)| 9 ) ~n"< 9) (30) 
that is expected to present a sort of phase transition [25[ . 

A. Continuous Time Random Walk approach 

For completeness, we briefly recall the standard th eory 
of Continuous Time Random Walks |U E 0, 0, Hf . 
Generally speaking, a CTRW is a stochastic model in 
which steps of a simple random walk take place at times 



tj, following some waiting time distribution. Mathemat- 
ically, it is asserted that a CTRW is a (non-Markovian) 
process subordinated to a random walk under the oper- 
ational time defined by the process fjij . 

A CTRW is completely characterized by the quan- 
tity ij)(r,T), the probability density function to move 
a distance r during a time interval r in a single mo- 
tion event; the dependence upon r and r can be ei- 
ther decoupled (i.e. ip(r,r) = x( r )p( r )) or coupled (e.g. 
ip(r,r) = x(r\r)p(T)). 

The object we are interested in, is the probability den- 
sity function P(x, t) of being in x at time t; indeed it 
allows us to obtain the full spectrum of transport mo- 
ments, through the formula 



(x(t)") = (i)«L- 



d q ~ 



(31) 



where is the inverse Laplace transform and P de- 
notes the Fourier-Laplace transform, being k the Fourier 
variable and u the Laplace variable. 
Let's introduce cf>(x,t), the probability density function 
of passing through (x,t), even without stopping at x, in 
a single motion event 

/>OC poo 

4>(x,t) =P(x\t) / dr drip(r,T). (32) 

Jt J\x\ 

P(x, t) is given by the sum of the probabilities of passing 
through (x, t), even without stopping at x, in one or more 
motion events 

/oo />£ 
dx' dri(j(x' > T)^(x-x',t-T) + ... 
-oo JO 

(33) 

By performing the convolution integrals, the Fourier- 
Laplace transform of this expression assumes the closed 
form 



P(k,u) 



1 — ip(k, u) 



(34) 



A special realization of CTRW is the so called velocity 
model 42]: a particle moves at a constant velocity for 
a given time, then stops and chooses a new direction 
and a new time of sojourn at random according to given 
probabilities. 

Our case belongs to this class, with velocities being ±1, 
and 



x (r|r) = ^8(\r\ - r) and p(r) 



(35) 



so that 



ip(r,r) ~ —8(\r\—T)T 9 and 



[X,t) -S(\ X \-t)t-" +1 

(36) 

where g = , being p(r) given by the waiting time 
distribution ipo,(n) of Eq. 
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By making use of the Tauberian theorems for the 
Laplace transform [4(| and by applying the CTRW for- 
malism [45[ we derive, through Eq. (f3Tj) and (|34|) the full 
spectrum of transport moments. 

The obtained spectrum of moments (more precisely, 
from the previous calculation it is possible to obtain only 
the even moments, and then to infer that a similar law 
drives also the behavior of the absolute value of odd mo- 
ments) is: 



(\z(n) - z(0)| 9 ) ~ 



(37) 



where the exponent v{q) has a picccwisc linear behavior 

7 + 3 



q/2 if q < 2a 
q — a if q > 2a 



7 



1 



in agreement with numerical results shown in Fig. ([5]). 
The transition at q — 2a in momenta spectrum of Eq. 



V(q) . 




FIG. 9: (color online) Spectrum of the transport moments 
for different value of the parameter 7. Lines correspond to 
theoretical predictions of Eq. (|38[1 . symbols correspond to 
numerical simulations: circles 7 = 3, diamonds 7 = 11/3, 
squares 7 = 5, triangles 7 = 7. 

(I38p is general in systems manifesting anomalous diffu- 
sion |2a |. 

As an outcome, we have that (anomalous) transport 
properties fully agree with the power laws we deduced 
for the waiting time distribution (Eq. (|26[) L 



VI. NOISE EFFECTS 

In order to better understand the link between cor- 
relations decay and time statistics (and to verify it, if 
not rigorously prove), we consider the effects of a small 
stochastic perturbation. The behavior, under the modi- 
fied dynamics, of the survival probability and of correla- 
tions decay may provide further informations about the 



interconnection between them. At the same time the dy- 
namical effects of a superimposed noise are interesting by 
themselves (see references (lt j l37 l l47l. l48j ). 
The general expectation is that small scale stochasticity 
blurs the behavior in the vicinity of the parabolic fixed 
point, enhancing the chaotic character of motion; one 
then expects a transition to an exponential decay of the 
survival probability and correlation functions: this intu- 
ition is corroborated by numerical experiments, reported 
in Fig. (fTX)|) . We perturb the system by introducing 



(38) C(n) 




2(H) 



FIG. 10: Correlation decay for noisy dynamics, for 7 = 10 
and various values of e (from top to bottom: e = 0.0, e = 
0.005, e = 0.010, e = 0.015, e = 0.020, e = 0.025, e = 0.030, 
e = 0.035, e = 0.040, e = 0.045, e = 0.045, e = 0.050). Each 
correlation function is computed by considering 3 • 10 9 initial 
conditions. 

a stochastic noise, adding at each iteration of the map 
a random vector of the type £ = (£1,^2) w hh £j i.i.d in 
(— e, e). The effects of the perturbation are expected to be 
dominant in the region of phase space (that will depend 
on the noise intensity e and on the stickiness parameter 7) 
where the deterministic step is small compared to noise. 
Basing on this assumption, we are able to give an ana- 
lytical estimate of the crossover time t c (defined as the 
characteristic time of the asymptotic exponential decay) , 
which is in good agreement with numerical simulations. 

Following [47J we divide the phase space into two com- 
plementary sections: one small region surrounding the 
fixed point, where the dynamics is dominated by stochas- 
tic diffusion, and its complementary, far from the fixed 
point, where the dynamics is dominated by the determin- 
istic chaotic motion. 

Now the main problem is a proper definition of the 
boundary of such a partition. The criterion suggested in 
|47| consists in defining (T ^deterministic) as the mean exit 
time from the region, say A, evaluated with the assump- 
tion that the dynamics is only due to the deterministic 
motion of the unperturbed map; then we define {t ra ndom) 
as the mean exit time evaluated as if the dynamics were 
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only due to diffusion. 

The border of the region is determined by the constraint: 



(Tdet erministic) — {^random} ■ 



(39) 



We restrict the analysis to the first quadrant and 
choose as region A the area defined by the survival prob- 
ability pn(fc), for some time k. In this way (T^ termimstic ) 
is given by 



(^deterministic) 



1 X 

— - -J2^n)-(n-k+l). (40) 



n>k 



Performing the calculation by substituting the probabil- 
ities obtained in the previous sections (Eq. (0 [5])), we 
get 



erministic/ 



k. 



(41) 



The calculation of (t ra ndom) is performed as follows: 
firstly we approximate pn(k) as in section III with the 
rectangular regions of Eq. ([24]) , say Xk-yk', these rectan- 
gular region can be exited along the x— direction or along 
the y— direction, independently (thanks to the particular 
form of our noise) , so that 



random! 



,.k,x v ,,k,y > 
V random' ' V raridom! 



(42) 



From the diffusion equation describing stochastic dynam- 
ics (see [13,lii]) we get 



random! 



(43) 



where z can be either x or y. Remembering that close to 
the origin the motion follows dynamics on the unstable 
manifold and from Eq. 



Xk 
Vk 



By substituting Eq. 



(^random) 



and finally, from Eq. ([39 



we get 

4 ~ k 

in Eq 



7-1 

, -2+1 

k -i- 1 . 



we obtain 



(fc-^T) 2 



(44) 



(45) 



(46) 



Writing the characteristic time of the exponential decay 
as a function of the noise strength e and the intermittency 
parameter 7, we derive an estimate for the crossover time 



tr 



-/3(7) 



27-2 

e 3 ^+ x . 



(47) 



Eventually we compare the behavior of ,8(7) with the 
numerical results, either for the survival probability and 
for the correlations decay. 



Numerical results are obtained in the following way: 
for each value of 7 we consider either the survival proba- 
bility or correlation functions for several values of e (see 
Fig. (fTU|): the exponential decay rates, that we consis- 
tently observe, are then fitted according to a power-law 
in e. 

The nice resemblance between the numerical simula- 
tions for the correlations decay and the survival proba- 
bility, together with a good agreement with the analytical 
result of Eq. (|4"T|) (see Fig. [TTj) strengthen our belief that 
the two distributions in the unperturbed case should be 
driven by the same exponent. 




FIG. 11: (color online) Comparison between theory and nu- 
merical results for the function /3(7) of Eq. (|47|l . Circles: nu- 
merical data for the survival probability; squares: numerical 
data for the correlations; dashed line: theoretical prediction. 



VII. CONCLUSIONS 

We considered a one-parameter family of area- 
preserving maps which generalizes intermittent behavior 
in two dimensions (while admitting a smooth invariant 
measure) . This is a paradigmatic example of weak chaos 
for hamiltonian maps even if the measure of the "regu- 
lar" portion of the phase space is zero (only a parabolic 
fixed point): the parameter 7 controls sticking of tra- 
jectories to the fixed point. By considering the motion 
along the unstable manifold close to the fixed point we 
are lead to estimates for power law decays of the sur- 
vival probability and waiting time distribution. Numer- 
ical computations of survival probabilities and residence 
time statistics show a very close agreement with pre- 
dicted exponents. Correlation functions are harder to 
deal with numerically, yet in this case results are close to 
analytic predictions. Our results are further supported 
either by considering transport properties for a lift of the 
map, and by studying dynamical effects induced by a 



9 



stochastic perturbation. 

While we think that a complete -quantitative- under- 
standing of weak chaos in more that one dimension is 
still in its infancy, our results shed new light on connect- 
ing local features (motion along the unstable manifold 
close to the fixed point) to global dynamical quantities, 
as mixing speed, transport properties, and response to 
stochastic perturbations. 
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